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Abstract 

Pick an arbitrary timelike geodesic in an arbitrary spacetime. We demon- 
strate that there is a particular limiting process, an "ultra-local limit", in 
which the immediate neighborhood of the timelike geodesic can be "blown 
up" to yield a general (typically non-diagonal) Bianchi type I spacetime. This 
process shares some (but definitely not all) of the features of the Penrose limit, 
whereby the immediate neighborhood of an arbitrary null geodesic is "blown 
up" to yield a pp-wave as a limit. 

Keywords: ultra- local limit, Bianchi type I spacetime, homogeneous anisotropic 
cosmology, Penrose limit, pp-wave. 

File: Bianchi-l-as-a-limit-27aug2010.tex 



1 



Contents 



1 Introduction |3j 

2 Ultra-local limit 

3 General Bianchi type I spacetime 5 

3.1 Ricci tensor and scalar 5 

3.2 Einstein tensor 6 

3.3 Mixed components § 

3.4 Unimodular decomposition 7 

3.5 Summary |8 

4 Stress-energy in the ultra-local limit 8 

4.1 Scalar field 8 

4.2 Electromagnetic field 19 

4.3 Perfect fluid ttO 

4.4 Anisotropic fluid 10 

4.5 Stress-energy conservation 11 

4.6 Summary Il2 

5 Einstein equations \12 

6 Conditions for a [block] -diagonal metric \13 

7 Timelike Raychaudhuri equation ll5 

8 Discussion 16 



2 



1 Introduction 



We shall demonstrate below that there is a particular limiting process, an "ultra-local 
limit", in which the immediate neighborhood of any arbitrary timelike geodesic of 
any arbitrary spacetime can be "blown up" to yield a general (typically non-diagonal) 
Bianchi type I spacetime. This ultra-local limiting process is inspired by the well- 
known "Penrose limit" [U EJ El HJ El E] , whereby the immediate neighborhood of an 
arbitrary null geodesic in an arbitrary spacetime is "blown up" to yield a pp-wave. 

While the ultra-local limit we construct below and the Penrose limit have some 
key features in common, there are also some significant differences — in particu- 
lar the "hereditary property" (see particularly reference [3]) fails for the ultra-local 
limit, which does significantly constrain its usefulness. On the other hand, given 
the enormous amount of attention that has been paid to Bianchi type I space- 
times for purely cosmological reasons, (for a selection of classic and recent sources 
see [51 El [71 El El [lOl HH [121 H31 H]), there is a lot of material available to help 
us understand the limiting Bianchi type I spacetimes. One particular caveat that 
should be borne in mind is that in a cosmological setting almost all authors rapidly 
restrict attention to Bianchi type I spacetimes with a diagonal metric — this is not 
appropriate when dealing with the ultra- local limit, and typically we must retain full 
generality and consider non-diagonal Bianchi type I spacetimes. 

Note further that a closely related limit has sometimes been used in developing 
mini-superspace models for canonical quantum gravity [151 HH], while a related ap- 
proximation underlies the so-called "velocity-dominated" cosmologies [HI HE] and 
the BKL singularities [19]. The focus herein is rather different, and we shall be more 
interested in the classical properties of the Bianchi type I spacetimes that emerge 
from this ultra- local limit. 

2 Ultra-local limit 

To define the ultra-local limit pick an arbitrary timelike geodesic 7 and, in some time- 
like tube surrounding the chosen curve, choose a specific set of adapted coordinates 
coordinates {t,x l ) such that: 

• The curve in question lies at x % = 0. 

• At x l = the coordinate t is just proper time along the curve. 

• The coordinate system is synchronous. 
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This coordinate choice is enough to set 

ds 2 = -c 2 dt 2 + g i:j (t, x) dxMx j , (1) 

where we have explicitly introduced the speed of light. Now implement the ultra-local 
limit in stages: 

1. Make the coordinate transformation x l — > ex 1 . Then 

ds 2 -c 2 dt 2 + e 2 g ij (t, ex) dx { dx j . (2) 

2. Change the speed of light c — > ec. So the net effect is 

ds 2 -e 2 c 2 dt 2 + e 2 gij (t, ex) dx*dx j . (3) 

3. Perform a conformal transformation ds 2 — > e~ 2 ds 2 . So overall 

ds 2 ->• -c 2 dt 2 + gij(t, ex) dx i dx j . (4) 

4. Take the limit e — > 0. Then the net effect is 

ds 2 -> -c 2 dt 2 + g tj (t) dx l dx j . (5) 

This definition is clearly "ultra-local" — one has effectively set c — > so that nothing 
propagates, so that every point in "space" becomes disconnected from every other 
point; and each individual point in space is then "blown up" by the conformal trans- 
formation to yield an independent universe. This closely parallels the Penrose limit 
construction [TJ [2j [3j H] , in fact one can argue that it is as close as one can possibly 
get to a "timelike Penrose limit" — though one should certainly not expect all of 
the features of the usual "lightlike Penrose limit" to survive in this timelike case. 

Note that the the 3-dimensional spatial slices are examples of the simplest (type I) 
element of Bianchi's classification of translation invariant 3-dimensional Riemannian 
geometries (20], so the output from this ultra- local limit is is exactly the definition 
of a general Bianchi type I spacetime [HI El d [HI 01 HQl Ell E2l H31 [H] — the general 
definition before you make any assumptions about diagonalizability of the spatial 
metric gij(t). Note further that if one starts with a congruence of timelike geodesies 
in the original spacetime, then the output of this ultra-local limit is a collection of 
general Bianchi type I spacetimes, one Bianchi spacetime being attached to each 
timelike geodesic. 
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3 General Bianchi type I spacetime 



Now that the speed of light c has done its job, let us quietly adopt units where c = 1, 
so that we are considering 

ds 2 = -dt 2 + gij (t) dsW. (6) 
Geometrically this is a rather simple spacetime. Define the extrinsic curvature 



= 4%, (7) 
3 2 dt 

then in particular 



HK) = --L (8) 

V93 dt 



where the notion of "trace" implicitly involves appropriate factors of the 3-metric 
and its inverse. Specializing the usual ADM decomposition, the key results for the 
Riemann tensor are 

Rijki = K ik Kji — KuK jk ; (9) 

and 

Roan = ^T + K imK m j = ^ + (K%; (10) 

with all other components being zero. For an explicit computation, see (and appro- 
priately modify) the discussion of the ADM decomposition in Misner, Thorne, and 
Wheeler [2T]. Note carefully that this ultra-local limit does not share the hereditary 
properties of the ordinary Penrose limit. In particular, the ultra- local limit of a 
Ricci flat spacetime is not necessarily Ricci flat — from the ADM view this is with 
hindsight obvious since during the limiting process the extrinsic curvatures and the 
intrinsic (spatial) curvatures scale in a different manner. 

3.1 Ricci tensor and scalar 

A brief computation yields 



while 



Therefore 



^ = ^- tr( A-), (11) 
at 



dK ■ 

Rii = -^-2K m K m J +tr(K)K lj . (12) 



R = - 2 ^^- + [tT(K)] 2 + tT(K 2 ). (13) 
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3.2 Einstein tensor 



The Einstein tensor can also be evaluated from the above, or can be read off from 
Misner, Thorne, and Wheeler [21]. See (21.162 a,b,c) p 552. The time-time compo- 
nent is easy 

Goo = i{[trW] 2 -tr(K 2 )}. (14) 
In contrast the space-space components are relatively messy 
Gij = ~ [K i:i - QiMK)] + - - \g t3 {MK)f + tr(K 2 )} . 



(15) 



3.3 Mixed components 



It is often advantageous to work with mixed components i?V Working from the 
above, (or suitably adapting results from the key article on BKL spacetimes |19j). 
one has: 

tf o = _dtgO +tr(x>) . (16) 

R'o = 0; (17) 
R°i = 0; (18) 

Ri '—itM^ Ki ' ] - (19) 



Also note that this implies (summing only over the space indices) 



*' = ^fS^ (20) 



For the Einstein tensor the interesting piece is 

1 d r _ _ . , „ fdtr(A') [tr(A')] 2 tr(A' 2 



dt 



(21) 
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3.4 Unimodular decomposition 

It is very useful to decompose 

9ij(t) = a{tf gij(t), 



(22) 



where det(^) = 1. Here a(t) can be viewed as an overall scale factor (similar to that 
occurring in FLRW cosmologies), while g^ describes the "shape" of space. Then 



and so 



— = 2aa 9ij + a —. 



-ad gij + a 2 Ki 



(23) 



(24) 



with = because of the unit determinant condition. Now define K l j = g lk Kkj 

so that 



Thus 



Then 



and 



a 2 



tr(K) = -3 -; tr(K 2 ) = 3 — + tr(K 



R° = 3- + tr(k 2 ); 
a 



R l i — — — 
dt 



K l 



a a z 



(25) 
(26) 
(27) 
(28) 



Summing over the space components only we have (in agreement with the previous 
subsection, see equation (120]) ) 



, a a 2 
3<! - + 2 — 

a er 



1 d 2 (a 3 ) 
a 3 dt 2 



(29) 



For the Einstein tensor the interesting piece is 



G\ = - — 
3 dt 



ly 3 



— 3-K\ - 



1 fa 



2 a 



+ 2- 



d 2 



(30) 
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3.5 Summary 



The spacetime curvature of the general Bianchi type I spacetime can be evaluated as 
simple algebraic combinations of the 3x3 matrices g^, g^, and g\j. If we split the 
geometry into an overall scale factor a(t) and a shape g^, then the curvature can be 
evaluated as simple algebraic combination of o, ci, d and the 3x3 matrices gij } gij-, 
and gij. 

4 Stress-energy in the ultra-local limit 

What happens to the stress-energy in the ultra-local limit we are interested in? Based 
on physical intuition one expects that for the energy- momentum flux Tq{ — > 0, but 
what about Too and Ty? (Remember that for the light-like Penrose limit T ab — > T uu 
with all other components vanishing [TJ |2j [3j H] . ) It is convenient to suppress explicit 
occurrences of c and to write the ultra-local limit in a slightly different form: 

1) Make the replacement 

g ab (t, x) -» g e ab (t, x) = g ab (t, ex). (31) 

2) Furthermore, for any generic field \I/(t, x) make the replacement 

(t, x) * e (t, x) = (t, ex). (32) 

3) Then consider the limit as e — > 0. 

4.1 Scalar field 

Consider a scalar field <f)(t, x) with stress-energy tensor 

T ab = 0, a 0, fe - l -g ab {(V0) 2 + V(<f>)} . (33) 
Now take the ultra- local limit. Note 

^oo = ~ {^e + gi^itedrfe + v } = i {0 2 + JjtdiWrf + v}^\ {0 2 (t, o) + 7 } , 

(34) 

whence ultimately 

^o {0 2 M) + ^} . (35) 
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Similarly 

= \ {kd4e} -+ \ {e&tf} -+ 0. (36) 

Furthermore 

Tt 3 = - \g\i { (V0) 2 + V} -+ e 2 ^-0 - I {(-0 2 + e 2 (90) 2 ) + \/} (37) 

whence ultimately 

T^{<p 2 (t,0)-V}g ir (38) 

The key point is that the structure of the stress-energy tensor guarantees that in the 
ultra-local limit 

T 00 ->■ p\ T 0i ->■ 0; Ty ->■ p ^. (39) 

This is a co-moving perfect fluid. Such simple behaviour is common, by no means 
universal. 

4.2 Electromagnetic field 

For the electromagnetic field the most natural form of the ultra-local limit is 

A a (t, x) ->■ A*(t, x) = ($(t, ex); A(t, exfj , (40) 

in which case 

E -> -A(t,0); B^O. (41) 
So for the stress-energy tensor 

Too^^j^M}; (42) 
T 0l 0; (43) 
Tij ^ - l -[g kl A k A^ 9ij (44) 

In terms of the electric field E, and the naturally defined inner product g(E,E), it 
is best to summarize this as: 

Too ^ <?(£,£); (45) 
T 0i -> 0; (46) 
T ij ^E i E j - l -g{E,E)g ij . (47) 



Note that in general the 3x3 matrix of space-space components of the stress-energy 
will not be diagonal. 

Note further that in reference [9] a somewhat different limit is taken — because 
those authors are interested in different physics. This is done because if one is 
interpreting the Bianchi type I spacetime as a cosmology, as an approximation to 
the large-scale structure of our own physical universe, then for physical reasons one 
might wish to allow B ^ 0, while expecting that E — > 0. Nevertheless, even in 
this very different context, one key point remains: In general the 3x3 matrix of 
space-space components of the stress-energy need not and will not be diagonal. 

4.3 Perfect fluid 

Consider a generic perfect fluid in a generic spacetime, 

T ab = (p + p)V a V h + pg a \ (48) 

and ask what the ultra-local limit might be? Working from the geometrical side we 
have G° l — > in the naturally adapted coordinate system. Thus via the Einstein 
equations it follows that we must also have T 0t — > 0, which in turn implies V 1 — > 0. 
Therefore 

Too ->■ P! T 0i -)■ 0; Tij -)■ p g^. (49) 

That is, the generic perfect fluid reduces in the ultra-local limit to a co-moving 
perfect fluid. 

4.4 Anisotropic fluid 

Consider a generic anisotropic fluid 

T ab = (p + p)V a V b + pg ab + Q ab , (50) 

and ask what the ultra-local limit might be? (Here is the anisotropic stress tensor, 
which is always 4-traceless and 4-orthogonal to the 4- velocity V.) However we choose 
to define this limit, since we know from the geometrical side that G° l — > 0, we must 
via the Einstein equations have have T 0t — > 0, so that the stress-energy tensor is 
(3+1) block diagonal. This in turn implies that we can without loss of generality 
redefine our variables so that V 1 — > 0, while V° — > 1, and simultaneously choose 
6 00 along with Oi ->- 0. Therefore 

Too -> P; T 0i -)• 0; -)• tt^ = pgij + 0^. (51) 
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4.5 Stress-energy conservation 

Consider the covariant conservation law 



V a T ah = —Lc> a (v^ T ab ) + Y h cd T cd = 0. (52) 



The only nontrivial component is 



! d t (^- 3 T 00 ) + T° cd T cd = 0, (53) 



which implies 



that is 



1 a/ ld 9^ij 



^ = -5-*^. (55) 

To see the connection with more usual FLRW spacetime note gij = a 2 and then 



3 d t {a 6 p) = -- [2aa 9ij - a 2 ) tt«, (56) 



which we can write as 



d t (a 3 p) = -a 2 a 9ij ^ - \a 6 ^-^ . (57) 
That is, (defining 7rjj = a 2 n^, so that 7r l j = if l j and 7r u = a' 2 7r lJ '), 

^!) = _ 3a v-i<» 3 *{*«-pn. (58) 

Finally, (defining 0^ = a 2 O -, so that 0^ = 0^ and Q ij = a~ 2 Q ij ), 



d (P° 3 ) = d ( Q3 ) _ I fl3 d&j Qi,' 



= -p^TLZ _ 1 a ^ ^ e tJ . (59) 
dt y dt 2 dt y ' 

The first 2 terms are the usual energy conservation law. The last term involves the 
trace-free anisotropic part of the (spatial) stress tensor, together with the trace-free 
contribution from dg^/dt. We can also write this as 

^ = -3(p + p)^ + £ y e*. (eo) 



11 



4.6 Summary 

In general the message is that in the ultra-local limit the stress-energy satisfies 

T o^p; T 0l ^0; T,, -> (61) 

and that one cannot say anything more than this without specifying the particular 
form of the matter content. Furthermore in this ultra-local limit there is a natural 
generalization of the energy conservation law normally applied to FLRW spacetimes, 
with an extra term coming from the interplay between anisotropies in the stress 
tensor and changes in the shape (not volume) of the spatial slices. 



5 Einstein equations 

In the ultra-local limit the Einstein equations, G a b = R a b — \Rg a b = 87rGArT a b reduce 
to three significant pieces of information — arising from the the 00 component, the 
spatial trace, and the anisotropic part of the spatial trace. It is most convenient to 
rewrite the Einstein equations as R a b = ^G^{T ab — ^Tg a b} and then adopt units so 
that Gn — > 1 in which case 

R 00 = 4tt( P + 3p); R\ = 12n(p - p); R tj - ^R k k 9ij = 8tt Q tj . (62) 

The first of these equations implies 

- = ^(p + 3p)-:W). (63 ) 

That is, as compared to spatially flat FLRW spacetimes, changes in the shape of the 
spatial slices can now contribute to deceleration. The second equation yields 

- + 2^- =47r(p-p). (64) 
a cr 

So this combination of terms is what one might expect anyway for any spatially flat 
FLRW cosmology. Finally 



d r 

dT 



a - 



3_ir i = 87reY (65) 



a 



That is, evolution of the shape of the spatial slices is driven by anisotropies in the 
stress tensor. This is the simplest and most straightforward decoupling of the Ein- 
stein equations we have managed to find. 
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6 Conditions for a [block] -diagonal metric 



A key issue in taking this ultra-local limit (and in Bianchi type I cosmologies gen- 
erally) is the question of when it is possible (or desirable) to restrict attention to 
diagonal spatial metrics (or more generally block diagonal metrics). Fortunately, 
within the context of Bianchi type I spacetimes it is possible to prove the following 
qualitative result: 

([block] diagonal stress-energy) <^ ([block] diagonal metric). (66) 



Proof (<=): A [block] -diagonal metric Qijif) implies Qijit) is [block]-diagonal, im- 
plies Kij(t) is [block]-diagonal, implies Kij(t) is [block]-diagonal, implies Rij(t) is 
[block]-diagonal, implies Gij(t) is [block]-diagonal implies TTij(t) is [block]-diagonal, 
implies 7r l j(t) is [block] -diagonal. ( Equivalent ly, this implies ®ij{t) is [block]-diagonal, 
which implies Q l j(t) is [block] -diagonal.) 

Proof (=>•): At time to there is no loss of generality in picking coordinates such 
that both [go]ij = and such that [go]ij is diagonal. Then without loss of generality 
there is a coordinate system such that [Ko]ij an d [-^oPj are diagonal. Our definition 
of [block]-diagonal stress-energy will be to assert that in this particular coordinate 
system T l j(t) is [block]-diagonal, whence G l j{t) is [block]-diagonal, and so R l j(t) is 
[block]-diagonal. But then 

^[V^l « {[block]-diagonal(*)} i j , (67) 

and so 

y/gilCj = ^[K^h + {[blockJ-diagonaKt)} 1 ,, (68) 
So by our initial assumption on the diagonal nature of [K ] l j we have 

ICj = {[block]-diagonal(t)}V (69) 

That is 

V* : g m g mj = {[block]-diagonal(*)}V (70) 

but then 

Vt : gij = g im {[blockj-diagonaKt)}" 1 ,, (71) 
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which can be formally integrated in terms of a time-ordered product 

\9(t)]ij = [tfoWexp [{[block]-diagonal (£)}]"%• (72) 

= 5i m exp [{[block]-diagonal(t)}] m j (73) 

= exp [{[block]-diagonal(t)}]^. (74) 

= {[block]-diagonal(t)}jj, (75) 

implying that [g(t)]ij remains [block] -diagonal for all time. 



Summary: A necessary and sufficient condition for a [block]-diagonal metric (in 
the context of Bianchi type I spacetimes) is that, in the coordinate system where 
[go\ij = $ij and [K ]ij is diagonal, we have 

n l j(t) is [block] -diagonal. (76) 

The metric and stress tensor then decompose into direct sums: 

9ij = ®A [h A }ij; TTij = ® A {p A [hA]ij}; = ®a {pa^aYj}- (77) 



Examples: In 3 space dimensions the three canonical examples of this behaviour 
(for a fully diagonal metric) are: 



1. 



Pefect fluid: 
That is: 



7T. 



P9ij- 



7T' 



(78) 
(79) 



p 
p 
p 

The special sub-case p = corresponds to dust, and the even more special 
sub-case p = p = corresponds to vacuum (Kasner solutions) — with all of 
these special cases (perfect fluid, dust, vacuum) being very well studied in a 
cosmological setting. 



2. Electromagnetic field with the electric field E being an eigenvector of g: 



This is equivalent to: 



7T, 



Pi gu o o 

P2 #22 o 
p 2 #33 



9ij 



gn o 

#22 








#33 



(80) 



(81) 
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That is: 





Pi 








ir* s = 


















P2 



3. General diagonal metric: 

7Ty = diag(pi,p 2 ,P3) 
This is equivalent to: 



7T 



'J 



Pig ii 
o 
o 



That is: 





P2 922 





7T 






P3 933 



(82) 



gij = diag(# u , g 22 , 533) 



<?n 
g 2 2 

£33 



(83) 



54) 



Pi 











P2 











P3 



(85) 



There is also the fully generic non-diagonalizable case corresponding to the stress 
tensor being a single irreducible 3x3 block. (And a partially diagonalizable case 
corresponding to one irreducible 2x2 block plus a singleton lxl block.) These are 
the only relevant cases for 3- dimensions. 



7 Timelike Raychaudhuri equation 

In a Bianchi type I spacetime we can simplify the standard (timelike) Raychaudhuri 
equation [22l [23] by noticing that we have a natural timelike congruence defined 
by it = d t , while we also have u = — (di)", or more prosaically u a = (1,0,0,0). 
For this timelike congruence the vorticity is zero, u = 0, as is the acceleration 
^j- = u h VbU a = 0. Therefore in terms of the expansion tensor ah = V( a w&), (not to 
be confused with the anisotropic stress Q a b), we have 

R ab u a u b = -e ab e ab + o 2 -^- (86) 

Separating out the unit determinant metric, = a 2 gij, and using the standard 
definitions of shear and expansion we find 

9 ab = aa g ab + a 2 K ab ; a ab = a 2 K ab ; a 2 = ^-ti(K 2 ). (87) 
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Furthermore K = 



—3d/ a, and hence 



a d8 a a 2 , . 

9 = 3-; — = 3- -3-2 88 
a at a a z 

So the Raychaudhuri equation is given by 

R ab u a u b = -tT(K 2 ) -3-. (89) 

(X 

Compare this to both the result for (spatially flat) FLRW cosmology, 

R ab u a u b = -3-, (90) 
a 

to which the Bianchi type I spacetime reduces when K ab = 0, and to the result stated 
by Collins and Ellis [B\ for perfect fluid Bianchi cosmologies, 

I(p + 3p) = -3^-2a 2 . (91) 

Attempting to put the Raychaudhuri scalar into a "nice" form, the best possible 
forms seem to be 

R ab u a u b =^(p + 3p)+tr(k 2 ), (92) 

or alternatively 

Ra b u a u h =-(p + ?,p)-l-. (93) 
4 2 a 



8 Discussion 



Overall the main message to take from this article is that the Bianchi type I space- 
times are interesting not only for their traditional use in cosmology, where they are 
the simplest example of homogeneous but anisotropic cosmologies, but also in the 
much wider context of the ultra-local limit described in this article. In a very precise 
manner described above, the Bianchi type I spacetimes are the result of a limiting 
process that can be applied to any timelike geodesic in any arbitrary spacetime. 
This ultra-local limit, although not sharing all important properties, seems to be the 
closest timelike analogue one can construct to the more usual Penrose limit (which 
is applied to lightlike geodesies to generate pp-wave spacetimes). 

In the current article we have first set up the ultra- local limit, and outlined its 
major features. One key item to note is that in the current context there is typically 
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no good reason to restrict oneself to diagonal Bianchi type I spacetimes — it is general 
(typically non-diagonal) Bianchi type I spacetimes that are of interest here. Modulo 
this technical issue, the existence of a vast cosmological literature on Bianchi type I 
spacetimes represents a significant resource which can be applied to the ultra-local 
limit of interest in this current article. 
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